YITP- 10-44 



O 

(N 

(N 
(N 

o 

u 

43 

Oh- 
i < 

o 

CO 

o3 



> 

oo 

O 

o 
o 



13 



Analytical model for CMB temperature angular power spectrum from cosmic 

(super-) strings 

Daisuke Yamauchi 1 !!] Keitaro Takahashi 2 Yuuiti Sendouda 3 Chul-Moon Yoo 1 ! and Misao Sasaki Q 
1 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
2 Department of Physics and Astrophysics, Nagoya University, Nagoya 494-8602, Japan 



APC, Universite Paris 7, 10, 



rue Alice Domon et Leonie Duquet, 
(Dated: June 23, 2010) 



75205 Paris cedex 13, France 



We present a new analytical method to calculate the small angle CMB temperature angular power 
spectrum due to cosmic (super-)string segments. In particular, using our method, we clarify the 
dependence on the intercommuting probability P. We find that the power spectrum is dominated 
by Poisson-distributed string segments. The power spectrum for a general value of P has a plateau 
on large angular scales and shows a power-law decrease on small angular scales. The resulting 
spectrum in the case of conventional cosmic strings is in very good agreement with the numerical 
result obtained by Fraisse et al.. Then we estimate the upper bound on the dimensionless tension 
of the string for various values of P by assuming that the fraction of the CMB power spectrum due 
to cosmic (super-)strings is less than ten percents at various angular scales up to I = 2000. We find 
that the amplitude of the spectrum increases as the intercommuting probability. As a consequence, 
strings with smaller intercommuting probabilities are found to be more tightly constrained. 



PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Cosmic strings are line-like topological defects formed 
in the early universe through spontaneous symmetry 
breaking in a wide range of inflationary models [l| . Since 
the string tension /i is directly related to the symmetry 
breaking energy scale, observational verification of the 
existence of cosmic strings will have profound implica- 
tions to unified theories. Theoretically, recent develop- 
ments in string cosmology suggest that inflation may be 
due to motions of branes in higher dimensions and vari- 
ous new types of strings, called cosmic superstrings, may 
be formed at the end of inflation One of the dif- 

ferences between cosmic superstrings and conventional 
field-theoretic strings is the value of the intercommuting 
probability P. It can be significantly smaller than unity 
for cosmic superstrings , while normally it is unity for 
field-theoretic strings @- 

It is known that a moving cosmic (super-)string in- 
duces a discontinuity in the gravitational potential, hence 
a discontinuity in the cosmic microwave background 
(CMB) temperature map. This is called the Gott-Kaiser- 
Stebbin (GKS) effect If photons are scattered by a 

number of moving string segments, the observed temper- 
ature fluctuations appear as a superposition of the dis- 
continuities. The imprint of cosmic strings on CMB has 
been widely studied. Using observed CMB anisotropy 
data, an upper bound on the dimensionless energy scale 
G/i was discussed The value varies but it is in 
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the range from 10 -7 to 10~ 6 (for P = 1). Although cos- 
mic strings were excluded as a dominant source of the 
observed large angle anisotropy, a signal due to cosmic 
strings could still be observed on small angular scales 
with future arcminutcs experiments such as South Pole 
Telescope [l4[ and Atacama Cosmology Telescope [lj|. 
A better theoretical understanding of the temperature 
fluctuations due to cosmic (super-)strings may help us 
distinguish them from other secondary effects and may 
enhance the observability of cosmic strings in such ex- 
periments. 

Recent numerical simulations [l6l - fl8| show that the 
small scale CMB temperature angular power spectrum 
due to cosmic strings with P = 1 behaves as a power 
law. One of our purposes of this paper is to derive this 
power-law behavior analytically and to extend it to the 
case of cosmic strings with P < 1 . 

In (l9l . |20| , we computed the one-point probability dis- 
tribution function (pdf) with a simple model of long 
curved string segments and kinks. It was found that 
the one-point pdf is dominated by a Gaussian component 
due to frequent scatterings by long straight segments, to- 
gether with small non-Gaussian tails due to close encoun- 
ters with kinks and a small asymmetry, namely skewncss, 
from the correlation between curvatures and velocities of 
string segments. Therefore, as far as the power spectrum 
is concerned, it is sufficient to consider only the contri- 
bution from long straight string segments. 

In this paper, we present a new analytical method to 
calculate the CMB temperature angular power spectrum 
due to cosmic (super-)string segments. Our formalism is 
similar to the halo formalism for the Sunyaev-Zcl'dovich 
effect [2ll - [23l | . We adopt a simple model of string network 
for general values of the intercommuting probability P. 
We consider long straight string segments which are lo- 
cated randomly between the last scattering surface (LSS) 
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and the present time consistently with the string network 
model. We find that the angular power spectrum is dom- 
inated by Poisson-distributed segments. Then we find it 
is possible to derive the dependence of the CMB spec- 
trum on the intercommuting probability P explicitly. 

This paper is organized as follows. In section lU we 
give basic equations governing a string network which 
incorporate the intercommuting probability P. In sec- 
tion IIII1 we give an explicit form of the GKS effect and 
perform the Fourier transformation analytically. In sec- 
tion IIVI we introduce a formalism called 'the segment 
formalism' and derive an analytical formula for the an- 
gular power spectrum due to cosmic (super-)strings. In 
section|V] we calculate an upper bound on the dimension- 
less tension G[i under appropriate assumptions. Finally, 
we summarize our results in IVII 



II. STRING NETWORK MODEL 

We first give basic equations governing string net- 
work incorporating intercommuting probability P. fol- 
lowing 0,H9]- A string worldsheet can be described by 
x n _ x n^ a a^ wnere x v anc j a a are t ne S p a cetime coordi- 
nates and the worldsheet coordinates, respectively. Let 
us consider string dynamics in a Friedmann-Lcmaitre- 
Robertson- Walker universe with the metric 



ds 2 =a 2 (r]) (-drf+dr 2 ) 



(1) 



We choose the temporal gauge : <r° = r\ , a 1 — a , r ■ 
r' = , where the bold letters denote 3- vectors on the 
comoving space and the dot and the prime denote the 
derivatives with respect to r\ and cr, respectively. 

In the velocity-dependent one-scale model, a string seg- 
ment has two properties, the length £ and the root-mean- 
squarc velocity v Tms . The typical length £ is defined by 
£ = \/ I 1 / Pscg where p scg is the total string energy density. 
In our treatment, we also take account of the energy loss 
due to loop formation. The characteristic time scale for 
loop formation is ~ £/(P f r ms) an d the energy loss can be 
described as ~ cPv Tms p seg /£ where we have introduced 
constant which represents the efficiency of loop 
formation. Assuming a(t) oc t@ with the physical time 
t = J a(rj)dri, the equations of motion for 7 and u rnls are 
given by HH HT 




p - \$~cPv n 



= C 1 - v 2 ms )H k{v T 



,2 

rms 1 



(2) 
(3) 



where k(v rms ) = (2y/2/n)(l - 8t&J/(l + 8z; r 6 ms ) M- 
Hereafter we assume a matter-dominated era, j3 = 2/3, 
and we use c ps 0.23 as the standard value (25j . 

It is known that a string network approaches the so- 
called scaling regime where the characteristic scale grows 
with the Hubble horizon size [111 [2?| ■ We assume that 
the scaling is already realized by the time of the last 



scattering surface (LSS) and this means that 7, u rms are 
constant in time. For small cP, we can solve Eqs. ©, 
(O approximately as 



7T 

37 



7 = l(P) 



! nV2 
3cP 



(4) 



III. 



TEMPERATURE FLUCTUATIONS DUE TO 
A STRING SEGMENT 



In order to discuss temperature fluctuations on small 
angular scales, we focus on a small patch of sky, and 
consider a straight string segment at the position r(r], a) 
where r\ and a are the time and position on the string 
worldsheet. We introduce a vector X{a) describing the 
comoving position of an observer relative to that of the 
string: 



X(p) = r ohs - r(cr, 7710(0-)) 



(5) 



where rj\ c {a) is the conformal time along the intersection 
of the observer's past light-cone and the string world- 
sheet, 7/ b s — ?/ic(c) = |-X"(er)|. On small scales, the tem- 
perature fluctuation due to a string segment in the direc- 
tion n, which is the unit vector along the line-of-sight, is 
given by [HQ 



AT 

— = -4G M 



da- 



X 



u 



X- 



(6) 



where X (a) = X(a) — (X(a) ■ n)n, and u is defined 
by 



u(a) = 



r — 



n ■ r 



1 



(7) 



and E denotes the intersection of the observer's past 
light-cone with the string worldsheet, along which the in- 
tegration is to be performed. We have adopted the small 
angle approximation and neglected higher order terms 
since their contributions are very small in general (20| . 

Hereafter, we assume that the string segment is exactly 
straight and uniformly moving, that is r' = const. ,r = 
const.. For an exactly straight and uniformly moving 
segment, Eq. (|6|) reduces to [2(1 [28| 



AT An H 



scg 



da 



(n x X- 1 ) 



da 



where we have introduced 



(8) 



(9) 



The position of the exactly straight and uniformly mov- 
ing segment can be written as 



aX ± (a) = d A "d + aa 



dX J 



da 



e. 



(10) 
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where <1a denotes the angular diameter distance from the 
observer, e oc dX 1 - /da = const, is the unit vector along 
the string, and # is the angular position vector relative 
to the middle point of the segment in a small patch of 
sky. 

Here, let us introduce a set of orthonormal basis vec- 
tors (si,S2,n), where the vectors (si,S2) span a patch 
of sky orthogonal to the line of sight. Then e and i? may 
be expressed as 



COS ifeSl + sin ip e S2 , 



i9 = d (cos ips\ + sin PS2) 



111) 



where the ranges of the coordinates are < # < oo and 
< tp < 2n in the flat sky approximation (which is valid 
under the small angle approximation). 

We also need explicit expressions of the components of 
the vectors r' and r, which determine the configuration of 
a string, where we adopt the gauge condition r' ■ r = 0. 
Thus the number of independent degrees of freedom is 
6—1 = 5. A convenient parametrization is (20| 



As mentioned above, since we focus on small angular 
scales, the flat sky approximation is valid. Therefore we 
can perform the Fourier transformation of the temper- 
ature fluctuation analytically. The Fourier transform is 
defined by (U S3 



(16) 



Then we obtain 



8iG/i|r|ow(0) . . / £ \ 

a * = /-, -9 ,9 — tan W sm ( T~ cos ft ) 1 ( 17 ) 

VI - r A l A \tco J 

where i = £(cos<^Si + sin<^s 2 ) and tpi = tpi — tp e - As 
apparent from the above expression, the result does not 
depend explicitly on tp e but only through the angle rel- 
ative to the string segment, tpi = pi — p e . This is the 
reason why it is unnecessary to express tp e explicitly in 
terms of the angular parameters 0. 



IV. SEGMENT FORMALISM 



r' = \r'\ sin 9 cos <frs i + sin 9 sin </>S2 + cos 9n , 

r = \r\ (— sin^ cosS cos0 — cos^sin0)si 
+ (— sin i\) cos 9 sin <f> + cos ijj cos 4>)s2 
+ sin ijj sin 9n , (12) 

where the ranges of the angular parameters are < 9 < 
7T , < (j) < 2?r, < -0 < 2tt. 

For notational simplicity, we introduce = {0,0,-0} 
to denote the set of the angular parameters. Then, 
Eq. (J5J) reduces to 



a sog (0) = cos tp sin 9 . 



(13) 



Note that the unit tangent vector e can be written as a 
function of the parameter 0, ip e = ip e (®). However, as 
we sec below, it turns out that we do not need an explicit 
expression for ip e . 

In our calculation, we consider only a segment of a 
long straight string with length ~ £ at each scatter- 
ing. Therefore we take the range of integration over a 
as a\dX^ / da\ \a\ < £. Then the temperature deviation 

is n ©,[ii 



AT \r\ 



x < arctan 



^co 1 + ^ cos <P 



i? sin tp 



arctan 



£ C o — d cos ip 



(14) 



where <p = tp — ip e and we have introduced the angular 
scale £ co corresponding to the correlation length of the 
segment, 



1 Hd A . 



(15) 



In order to compute the angular power spectrum of the 
temperature fluctuations due to cosmic (supcr-)strings, 
we use what we call the segment formalism, by adapting 
from the halo formalism for the Sunyaev-Zel'dovich ef- 
fect [2lTj23l | . Since the observed sky map of temperature 
fluctuations due to segments appears as a superposition 
of those due to each segment, the Fourier transform of 
the total temperature fluctuations, a^ ot , can be decom- 
posed into each contribution of each string segment. In 
our treatment, we first introduce a segment index "i" to 
denote the contribution from each segment between LSS 
and the present. Then we have 



/v 



if t ({& i ,z i }) = J2^(®i,z i ) 



1=1 



N 

t=i 



(i) 



(18) 



where 0^ and Zj are the segment configuration parame- 
ters and the redshift, respectively, of the i-th string seg- 
ment. N is the total number of the string segments. 

Asuming the statistical isotropy of the CMB, the an- 
gular power spectrum can be written as 



Ct = 



dipt/ 
2tt \ 
dipt 
~2^ 

lseg 



t ({0. ^ ,z ^ })ar({0„^}) , 



W| 2 



dpi 
"2V 



E« 



C, 



2scg 



(19) 



where the integral over pg is the large £ approximation of 
the sum over the azimuthal eigenvalues m (—£ < m < £), 
and (• • ■ ) denotes the ensemble average. The ensemble 
average can be calculated by averaging over the parame- 
ter space, 



n 



1 l A dV 



d® t 



dn 
d®~ 



(20) 
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FIG. 1: The contribution to Ce from each logarithmic interval 
of 1 + z, £(£ + l)D e (z)/(2ir), in the case of P = 1 in units of 
(Gfj,) 2 . Here it is shown as a function of £ for various values 
of z. From top to bottom, z = zlss, 500, 100, 10, 1 and 0.1. 



where [dV j dz)dz is the differential comoving volume el- 
ement at rcdshift z, (dn/d@) ■ d® is the comoving num- 
ber density of string segments with the parameters in the 
range [©, 0+d0] . Note that the total number of the seg- 
ments can be rewritten as N = J dz(dV/dz) J(dn/d®) ■ 
d®. We assume uniform distributions P(cos8) = 1/2, 
P(4>) = P(ip) = 1/2-7T of the parameters. Then, the num- 
ber density of the segments can be estimated as 



dn 
d®~ 



d&i 



3 3 d(coa9 i )d(f>idip i 



2(2tt) 2 



(21) 



where we have adopted the scaling ansatz and used 
Eq. ©. 

As mentioned in Sec. HI long straight string segments 
are assumed to be distributed randomly between LSS and 
the present consistently with the string network model. 
This implies there is no correlation between two seg- 
ments, (ai^ a'P) = for i =/= j. If we consider a more gen- 
eral string network, there may be some nonzero contri- 
bution from the segment-segment correlation. However, 
the segment-segment correlation at a rcdshift z would be 
dominated by the contribution from £ w £ co {z) and the 
smallest scale at which the segment- segment correlation 
could be important is determined by £ co at z = zlss, 



4o(z) < 4o(zlss) 



156 



023 



-1/2 



(22) 



where we have put zlss ~ 1100. Therefore, the angular 
power spectrum on small scales, £ > 4o(zlss)j will be 
dominated formally by the contribution of the sum of N 
statistically independent segments, even if the segment- 
segment correlation is taken into account. That is, 



Ct 



ZLSS d y 

dz—— 

dz 



dn 

dQ ._ g t (e t z), (23) 



with 



Gi{®,z) 



— \a £ (@,z)\ . 

Z7T 



(24) 



FIG. 2: The same as Fig. Q] but as a function of z for various 
values of £. From top to bottom, £ = 10 2 , 10 3 , 10 4 and 5x 10 4 . 



where the integral over ipi has been replaced by that over 
tf£ — > ipi = — tp e without loss of generality. 

Assuming that \r\ — v rms and that the universe is 
matter-dominated, we have £ co {z) w 27(-\/l + z — 1) and 
dV/dz = 16nH- 3 (l - (1 + z)- 1 / 2 ) 2 . Then we obtain an 
explicit form of Ci as 



1 + z 



(25) 



where 



D t (z) = (i + z)^ J d®--^g,(®,z) 



2^(8 7 ) 3 « r 2 ms (G/i) 



3(1 - v? ms )£ 



2n 



2 ^ 2 

tan top sin 



£co(z) 



cos tpi I . (26) 



Here Dg(z) is the redshift distribution of C*£, which tells 
us which redshift z contributes most for a given £, or 
which £ contributes most at a given z. In Figs. [1] and [21 
we plot Di(z) for P = 1. We see that the large z contri- 
bution dominates for a given £. This may be explained 
by the fact that the number of segments per unit red- 
shift bin becomes large as z increases. In other words, 
the contribution from small z is negligible, justifying the 
use of the small angle approximation in Eq. ([6]). 

The total CMB temperature angular power spectrum 
in the case of P = 1 is shown in Fig. [3l As seen from 
it, a typical amplitude of the power spectrum at £ = 10 
is [£(£+1)/2-k]C(£ = 10 3 ) m 14 (G^t) 2 -and it behaves as 
£- x for large £(> 4o(zlss)) (see also [MlHHl) while it 
has a plateau for small £(< 4o(zlss))- For comparison, 
we also plot power-law fit to previous numerical result: 



1) 



C t « 14(G^) 5 



( 



-0.89 



2tt v ^ V 1000 / 

for 400 < £ < 10 4 (Fraisse et al. Q] 
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FIG. 3: The angular power spectrum for P = 1 in units of 
(Gfj,) 2 . The red solid line is our result given by Eq. Q25p . The 
red dotted line shows the behavior oc The gray solid 

lines are power-law fit to previous numerical result by Fraisse 
et al. QJ]. 
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FIG. 5: The angular power spectrum in units of fj,K 2 . The 
red lines are for P = 1 and the blue lines are P — 10~ 3 . For 
both cases Gfi = 5 x 10" 7 , 10" 7 and 1(T 8 fr om top to bottom. 
For comparison, the primary spectrum is shown in gray. 




FIG. 4: The angular power spectrum given by Eq. (|25[) in 
units of (Gfj,) 2 . The curves are, from bottom to top, for P — 1 
(red), P = 10" 1 (orange), P = 10" 3 (green), and P = 10" 6 
(blue). 



As clearly seen, our result agrees very well with the 
numerical result by Fraisse et al. [l6j . This strongly 
supports the validity of our approach. We should note 
that the spectrum obtained here as well as the one 
obtained by Fraisse et al. [HI ignores fluctuations in- 
duced by strings on the last scattering surface (LSS). On 
small angular scales, the LSS contributions are damped 
and the GKS contributions dominate. Nevertheless the 
LSS contributions are not completely negligible even 
at i sa 4000, which seems to be the case in [l7l . Il8j |. 
Thus, one cannot really directly compare the spectra 
from Bevis et al. fl7l and Pogosian et al. [l8| using the 
code CMBACT |35( to those of our result and Fraisse 
et al. for I < 3000 0. In 0, the authors showed 
that the £~ 2 behavior for 1000 < £ < 3000 becomes 
much closer to £~ x when only the string sources after 
recombination are taken into account. Therefore, our 



as only the GKS component is considered. 1 

In order to investigate the dependence on the intercom- 
muting probability P, the angular power spectrum (|25p 
is computed for various P. The results are shown in 
Fig. [U We see that the overall amplitude of the spec- 
trum increases as P decreases. This is because of the 
factor 7 3 oc P~ 3 / 2 in the formula (|25|) . which describes 
the fact that the density of cosmic string segments is 
larger for smaller P. Also since £ co {zlss) °c P -1 ^ 2 , we 
see that the transition from the plateau to the power law 
occurs at larger £ for smaller P. These properties of the 
power spectrum may become a useful tool to distinguish 
the value of P in future experiments. 



V. CONSTRAINTS ON STRING TENSION 

Let us discuss possible constraints on the string tension 
from our result. We plot the angular power spectrum for 
various values of P and Gfj, in Fig. For comparison, we 
also plot the primary spectrum. An interesting observa- 
tion is that as P decreases the amplitude due to strings 
increases, hence the tension of strings with smaller P is 
more tightly constrained. 

It was pointed out in [10( that the CMB anisotropy 
spectrum is consistent with the presence of cosmic strings 
if the fraction of the power spectrum due to cosmic 
strings is about 10% or less at I = 10. In this paper 
we adopt a similar criterion and drive an upper bound 



result agrees not only with [16( but also with [17[ as far 



1 Our result does not agree well with Pogosian et al. |18( . They 
used CMBACT [35| and obtained a larger amplitude for 100 < 
£ < 3000 and an approximate i~ 1 ' 8 decay at high I. The differ- 
ence in the amplitude is probably due to different scaling param- 
eters such as 7 and « rma between theirs and ours. 
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FIG. 6: Upper bound on Gfi as a function of P at £ = 10 2 , 
5 x 10 2 , 10 3 and 2 x 10 3 , assuming that the fraction of the 
spectrum due to cosmic (super-)strings is less than 10%. 



on Gfj, as a function of P. Specifically, we consider the 
condition that the fraction of the power spectrum due 
to cosmic (super-)strings is less than 10% at I = 10 2 , 
5 x 10 2 , 10 3 and 2 x 10 3 . The result is shown in Fig. M 
As expected, the upper bound on G\i decreases as P de- 
creases, because the amplitude of the power spectrum in- 
creases. Also, we see that the constraint becomes severer 
for larger £ because the contribution from cosmic strings 
decays very slowly as £ increases unlike the case of the 
primordial anisotropy which shows exponential dumping. 
For example, the upper bound at £ = 10 2 is 2.1 x 10 6 for 
P = 1 and 3.1 x 1CT 8 for P = lfT 6 , while that at t = 10 3 
is 9.8 x lfT 7 for P = 1 and 6.0 x 10~ 9 for P = 10~ 6 . 



VI. SUMMARY 

In this paper, we presented a new analytical method 
to calculate the small angle CMB power spectrum due to 
cosmic (super-)string segments, and investigated the de- 
pendence of the power spectrum on the intercommuting 
probability P. 

We found that the angular power spectrum on small 
scales can be well approximated by the GKS effect due 
to Poisson-distributcd mutually independent segments. 
Then wc derived an analytical formula for the power 
spectrum valid for general values of P. The obtained 
power spectrum for P = 1, that is for conventional cos- 
mic strings, was found to agree very well with the numer- 
ical result obtained by Fraisse et al. [l6|. This strongly 
supports the validity of our approach, hence allows us to 
discuss the dependence on P with high confidence. 

The angular power spectrum is found to behave as £~ 1 
for large £(^> ^co(^lss)) and have a plateau for small 
^0$ ^co(^lss)), where ^co(^lss) is the angular scale cor- 
responding the correlation length at LSS. Since ^ C o(£lss) 
is proportional to P -1 / 2 in the scaling regime, the tran- 
sition from a plateau to the power-law behavior is found 



to occur at larger £ as P decreases. Wc should note, 
however, that the plateau region of the power spectrum 
may have additional modifications because the segment- 
segment correlatoin may not be negligible on scales £ < 
£ co (zlss)j though its effect on the spectrum is expected 
to be small if not negligible. 

Then using our result, wc discussed an upper bound 
on the dimensionless tension G\x as a function of P. We 
assumed that the fraction of the CMB spectrum due to 
cosmic (super-)strings is less than 10%, and derived an 
upper bound at I = 10 2 , 5 x 10 2 , 10 3 and 2 x 10 3 . We 
found that strings with small P are more tightly con- 
strained. This can be naturally explained by the fact 
that the amplitude of the spectrum increases as P de- 
creases because of the increase in the number density of 
strings. 

These properties of the power spectrum are distin- 
guishable features of cosmic superstrings that generally 
have a small intercommuting probability P. They may 
be used to detect cosmic superstrings in future experi- 
ments. 

Finally, wc comment on the bispectrum due to string 
segments. It is easy to see that in the present approach 
wc have ((a sog ) 2m+1 ) = and ((a sog ) 2m ) ^ for to = 
0, 1, • • ■ . This implies a vanishing bispectrum. To obtain 
a non-vanishing bispectrum it is necessary to take the 
correlation between the velocity and the curvature of a 
string segment into account [H, HH (see also 0, [2(| 
for discussion on the skewness of a one-point probability 
distribution function) . It is left for future work to include 
such correlations in our segment formalism, and calculate 
the non-Gaussianity of the spectrum (e.g. [13] )• Since 
the non-Gaussian features are expected to significantly 
depend on the intercomuting probability, they may be 
used to distinguish cosmic superstrings from conventional 
field theoretic cosmic strings. 
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